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Cosmic Ray Diffusion from the Galactic Spiral Arms, 
Iron Meteorites, and a possible climatic connection? 
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We construct a Galactic cosmic ray (CR) diffusion model while considering that CR sources 
reside predominantly in the Galactic spiral arms. We find that the CR flux (CRF) reaching the 
solar system should periodically increase each crossing of a Galactic spiral arm. We search for this 
signal in the CR exposure age record of Iron meteorites and confirm this prediction. We then check 
the hypothesis that climate, and in particular the temperature, is affected by the CRF to the extent 
that glaciations can be induced or completely hindered by possible climatic variations. We find that 
although the geological evidence for the occurrence of IAEs in the past Eon is not unequivocal, it 
appears to have a nontrivial correlation with the spiral arm crossings — agreeing in period and phase. 
Thus, a better timing study of glaciations could either confirm this result as an explanation to the 
occurrence of IAEs or refute a CRF climatic connection. 

PACS numbers: 98.35.Hj, 92.40.Cy, 92.70.Gt, 98.70.Sa 
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With the possible exception of extremely high ener- 
gies, cosmic rays (CRs) are believed to originate from 
supernova (SN) remnants (l], ^]. Moreover, most SNe 
in spiral galaxies like our own are those which originate 
from massive stars, thus, they predominantly reside in 
the spiral arms, where most massive stars are born and 
shortly thereafter explode as SNe ||. Indeed, high con- 
trasts in the non-thermal radio emission are observed be- 
tween the spiral arms and disks of external galaxies. As- 
suming equipartition between the CR energy density and 
the magnetic field, a CR density contrast can be inferred. 
It can have a lower limit of 5 in some cases pj . 

Thus, while the Sun is crossing the Galactic spiral 
arms, the CRF is expected to be higher. To estimate 
the CRF variation, we construct a simple diffusion model 
which considers that the CR sources reside in the Galac- 
tic spiral arms. We expand the basic CR diffusion models 
(e.g., ref. to include a source distribution located in 
the Galactic spiral arms. Namely, we replace a homoge- 
neous disk with an arm geometry as given by Taylor & 
Cordes M , and solve the time dependent diffusion prob- 
lem. To take into account the "Orion spur" pj , in which 
the Sun currently resides, we add an arm "segment" at 
our present location. Since the density of HII regions in 
this spur is roughly half of the density in the real nearby 
arms || , we assume it to have half the typical CR sources 
as the main arms. We integrate the CR sources assum- 
ing a diffusion coefficient of D = 10 28 cm 2 /sec, which is 
a typical value obtained in diffusion models for the CRs 

0, D . We also assume a halo half- width of 2kpc, which 
again is a typical value obtained in diffusion models 
but more importantly, we reproduce with it the 10 Be sur- 
vival fraction ||. Thus, the only free parameter in the 
model is the angular velocity Af2 = fi — tt p around 
the Galaxy of the solar system relative to the Spiral arm 
pattern speed, which is later adopted using observations. 
Results of the model are depicted in fig. 1. For the nomi- 
nal values chosen in our diffusion model and the particu- 



lar pattern speed which will soon be shown to fit various 
data, the expected CRF changes from about 25% of the 
current day CRF to about 135%. Moreover, the average 
CRF obtained in units of today's CRF is 76%. This is 
consistent with measurements showing that the average 
CRF over the period 150-700 Myr before presenWBP), 
was about 28% lower than the current day CRF Eg ]. 

Interestingly, the temporal behavior is both skewed 
and lagging after the spiral arm passages. The lag arises 
because the spiral arms are defined through the free elec- 
tron distribution. However the CRs are emitted from SNe 
which on average occur roughly 15 Myr after the aver- 
age ionizing photons are emitted. The skewness arises 
because it takes time for the CRs to diffuse after they 
are emitted. As a result, before the region of a given 
star reaches an arm, the CR density is low since no CRs 
were recently injected in that region and the sole flux is 
of CRs that succeed to diffuse to the region from large 
distances. After the region crosses the spiral arm, the CR 
density is larger since locally there was a recent injection 
of new CRs which only slowly disperse. This typically in- 
troduces a 10 Myr lag in the flux, totaling about 25 Myr 
with the SN delay. This lag is actually observed in the 
synchrotron emission from M51, which shows a peaked 
emission trailing the spiral arms y] . 

The spiral pattern speed of the Milky Way has not yet 
been reasonably determined through astronomical obser- 
vations. Nevertheless, a survey of the literature reveals 
that almost all observational determinations cluster ei- 
ther around Afl « 9 to 13 (km s -1 )/kpc|ll[] or around 
Afl k 2 to 5 (km s~ 1 )/kpc[Q. In fact, one analysis 
]l3| revealed that both AO = 5 or 11.5 (km s _1 )/kpc fit 
the data. However, if the spiral arms are a density wave 
|l4j| , as is commonly believed [|l5|, then the observations 
of the 4-arm spiral structure in HI outside the Galactic 
solar orbit EJI severely constrain the pattern speed to 
AJ7 > 9.1 ± 2.4 (km s _1 )/kpc, since the four arm den- 
sity wave spiral cannot extend beyond the outer 4 to 1 
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Lindblad resonance |L7j. We therefore expect the spiral 
pattern speed obtained to coincide with one of the two 
aforementioned ranges, with a strong theoretical argu- 
mentation favoring the first range. 

To validate the above prediction, that the CRF varied 
periodically, we require a direct "historic" record from 
which the actual time dependence of the CRF can be 
extracted. To find this record, we take a compilation of 
74 Iron meteorites which were 41 K/ 40 K exposure dated 
JlSl . CRF exposure dating (which measures the duration 
a given meteorite was exposed to CRs) assumes that the 
CRF history was constant, such that a linear change in 
the integrated flux corresponds to a linear change in age. 
However, if the CRF is variable, the apparent exposure 
age will be distorted. Long periods during which the 
CRF is low would correspond to slow increases in the ex- 
posure age. Consequently, Fe meteorites with real ages 
within this low CRF period would cluster together since 
they will not have significantly different integrated expo- 
sures. Periods with higher CRFs will have the opposite 
effect and spread apart the exposure ages of meteorites. 
To avoid real clustering in the data (due to one parent 
body generating many meteorites), we remove all occur- 
rences of Fe meteorites of the same classification that are 
separated by less than 100 Myr and replace them by the 
average. This leaves us with 42 meteorites. 

From inspection of fig. 1, it appears that the meteorites 
cluster with a period of 143 ± 10 Myr, or equivalently, 
|AQ| = 11.0 ± 0.8 (km s _1 )/kpc, which falls within the 
preferred range for the spiral arm pattern speed. If we 
fold the CR exposure ages over this period, we obtain 
the histogram in fig. 2. A K-S test yields a probability 
of 1.2% for generating this non- uniform signal from a 
uniform distribution. Moreover, fig. 2 also describes the 
prediction from the CR diffusion model. We see that the 
clustering is not in phase with the spiral arm crossing, but 
is with the correct phase and shape predicted by the CR 
model using the above pattern speed. A K-S test yields 
a 90% probability for generating it from the CR model 
distribution. Thus, we safely conclude that spiral arm 
passages modulate the CRF with a ~ 143 Myr period. 

In 1959, Ney |t9) suggested that the Galactic CR flux 
(CRF) reaching Earth could be affecting the climate since 
the CRF governs the ionization of the lower atmosphere, 
to which the climate may in principle be sensitive to. If 
this hypothesis is correct, we may be able to see a cor- 
relation between the observed long term CRF variability 
and the climate record on Earth. 

Interestingly, the CRF reaching Earth is also variable 
because of its interaction with the variable solar wind. 
Thus, solar activity variations will too have climatic ef- 
fects if the CRF affects the climate (e.g., Under the 
assumption that it does affect climate, we can estimate 
how large an effect can a possible CRF-temperature rela- 
tion be. This can be derived from the fact that the best 
fit to the global warming in the past 120 years is obtained 
if somewhat less then half is attributed to anthropogenic 
greenhouse gases and somewhat more than half to the in- 



creased activity of the sun pi, 22 . Thus, between about 
1940 and 1970, the global temperature, which decreased 
by 0.15°K, is best explained as — 0.2°K attributed to the 
reduced solar activity and +0.05°K to greenhouse gases 
pl| , [22] | . A global CRF climate effect is presumably more 
likey to arise from CRs that can reach the troposphere 
and equatorial latitudes. Thus, it is reasonable to as- 
sume that a possible effect would arise from CRs that 
have high rigidities (> 10-15 GeV/nucleon). We there- 
fore normalize the low geomagnetic data from Haleakala, 
Hawaii and Huancayo, Peru to the higher geomagnetic 
data of Cli max, Colorado |^3| that was measured over a 
longer period (e.g., |M|). We find that the -0.2°K cool- 
ing correlated with a 1.5% increase in the high rigidity 
CRF. Thus, changing the CRF by ±1% would correspond 
to a global change of =F0.13°K, on condition that CRs are 
indeed the link relating solar activity to the climate. 

For the nominal values chosen in our diffusion model, 
the expected CRF changes from about 25% of the current 
day CRF to about 135%. This corresponds to a tempera- 
ture change of about +10°K to — 5°K, relative to today's 
temperature. This range is sufficient to markably help or 
hinder Earth from entering an IAE. 

Extensive summaries of IAEs on Earth can be found 
in Crowell |^5| and Frakes et al. |^6| . Those of the past 
Eon are summarized in fig. 1. The nature of some of the 
IAEs is well understood while others are sketchy in detail. 
The main uncertainties are noted in fig. 1. For example, 
it is unclear to what extent can the milder mid-Mesozoic 
glaciations be placed on the same footing as other IAEs, 
nor is it clear to what extent can the period around 700 
Myr BP be called a warm period since glaciations were 
present, though probably not to the same extent as the 
periods before or after. Thus, Crowell (25) concludes that 
the evidence is insufficient to claim a periodicity. On the 
other hand, Williams J27J claimed that a periodicity may 
be present. This was significantly elaborated upon by 
Frakes et al. [^6| . 

Comparison between the CRF and the glaciations in 
the past 1 Gyr shows a compelling correlation (fig. 1). To 
quantify this correlation, we perform a x 2 analysis. To 
be conservative, we do so with the Crowell data which 
is less regular. Also, we do not consider the possible 
IAE around 900 Myr, though it does correlate with a 
spiral arm crossing. For a given pattern speed, we pre- 
dict the location of the spiral arms using the model. We 
find that a minimum is obtained for Afi = 10.9 ± 0.25 
(km s _1 )/kpc, with Xmin = 1-1 P er degree of freedom 
(of which there are 5=6-1). We also repeat the analy- 
sis when we neglect the lag and again when we assume 
that the spiral arms are separated by 90° (as opposed 
to the somewhat asymmetric location obtained by Tay- 
lor and Cordes d)). Both assumptions degrade the fit 
(Xmin = 2-9 with no lag, and Xmin ~ ^.1 with a symmet- 
ric arm location). Thus, the latter analysis assures that 
IAEs are more likely to be related to the spiral arms and 
not a more periodic phenomena, while the former helps 
assure that the CRs are more likely to be the cause, since 
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FIG. 1: The past Eon. Panel A describes past Galactic spi- 
ral arms crossings assuming AS7 = 10.9 (km s _1 )/kpc. Panel 
B describes the CRF reaching the solar system using the CR 
diffusion model. Note that the CRF lags behind spiral arm 
crossings. This is portrayed by the hatched regions, which 
qualitatively show the predicted occurrences of IAEs if the 
CRF required to trigger them is the average flux. Arrows 
mark the middle of the spiral crossing and the expected mid- 
glaciation point. Panel C qualitatively describes the geologi- 
cally recorded IAEs — its top half, as summarized by Crowell 
[ p5[ , who concluded that the current evidence is insufficient 
to support a claim of periodicity, while the bottom half, as 
summarized by Frakes et al. J2(| who claim a periodicity ex- 
ists. By fine tuning the observed pattern speed of the arms to 
best fit the IAEs, an intriguing correlation appears between 
the IAEs and their prediction. Note that the correlation need 
not be absolute since additional factors may affect the cli- 
mate. Other factors that should be considered and noted in 
the graph are: (1) The mid- Mesozoic glaciations are signifi- 
cantly less extensive as others. (2) It is unclear to what extent 
was the period around 700 Myr BP warmer than the IAEs 
before or after. (3) The first IAE of the Neo-Proterozoic (if 
indeed distinct) is very uncertain. (4) Since Norma's crossing 
is an extrapolation from smaller galactic radii, its location is 
uncertain. If the arm's structure at smaller radii is indeed 
different its preferred location will lag by about 20 Myr. 
Panel D is a 1-2-1 averaged histogram of the 41 K/ 40 K ex- 
posure ages of Fe Meteorites, which are predicted to cluster 
around the CRF minima. The cluster-IAE correlation further 
suggests an extra-terrestrial trigger for the glaciations. 



they are predicted (and observed) to be lagged. 

The previous analysis shows that to within the limi- 
tation of the uncertainties in the IAEs, the predictions 
of the CR diffusion model and the actual occurrences of 
IAE are consistent. To understand the significance of 
the result, we should also ask the question what is the 
probability that a random distribution of IAEs could gen- 
erate a x" 2 result which is as small as previously obtained. 
To do so, glaciation epochs where randomly chosen. To 
mimic the effect that nearby glaciations might appear as 
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FIG. 2: A folded 1-2-1 averaged histogram (over the IAE pe- 
riodicity) of the Fe Meteorites' exposure ages. A statistically 
significant signal is obtained for a period of 143 ±10 Myr. The 
dashed line is the CR diffusion model prediction. We find that 
the predicted phase of the cluster peaks, agrees with the ac- 
tual exposure age clustering. The clustering peaks 100 Myr 
after the spiral crossing (dotted lines), implying that it is un- 
likely to be a real age clustering peak, which should be related 
directly to the spiral crossing (e.g., by Oort cloud comets be- 
ing perturbed during arm crossing, that subsequently break 
asteroids more frequently, and that somehow affect the ter- 
restrial climate). 



one epoch, we bunch together glaciations that are sepa- 
rated by less than 60 Myrs (which is roughly the small- 
est separation between observed glaciations epochs). The 
fraction of random configurations that surpass the x 2 ob- 
tained for the best fit found before is of order 0.1% for any 
pattern speed. (If glaciations are not bunched, the frac- 
tion is about 100 times smaller, while it is about 5 times 
larger if the criterion for bunching is a separation of 100 
Myrs or less). The fraction becomes roughly 6 x 10 -5 (or 
a 4-er fluctuation), to coincidcntally fit the actual period 
seen in the Iron meteorites. 

Last, before 1 Gyr BP, there are no indications for any 
IAEs, except for periods around 2 - 2.5 Gyr BP (Huro- 
nian) and 3 Gyr BP (late Archean) pB| . This too has a 
good explanation within the picture presented. Differ- 
ent estimates to the Star formation rate (SFR) in the 
Milky Way (and therefore also to the CR production) 
point to a peak around 300 Myr BP, a significant dip 
between 1 and 2 Gyr BP (about a third of today's SFR) 
and a most significant peak at 2-3 Gyr BP (about twice 
as today's SFR)|2|, ||. This would imply that at 300 
Myr BP, a more prominent IAE should have occurred — 
explaining the large extent of the Carboniferous-Permian 
IAE. Between 1 and 2 Gyr BP, there should have been 
no glaciations and indeed none were seen. Last, IAEs 
should have also occurred 2 to 3 Gyr BP, which explains 
the Huronian and late- Archean IAEs. 

To conclude, by considering that most CR sources re- 
side in the Galactic spiral arms, we predict a variable 
CRF. A record of this signal was indeed found in Iron me- 
teorites, and it nicely agrees with the observations of the 
Galactic spiral arm pattern speed. Next, if the apparent 
solar activity climate correlation is real and arises from 
modulation of the galactic CRF reaching Earth, then typ- 
ical variations of up to O(10°K) could be expected from 
the variable CRF. Each spiral arm crossing, the aver- 
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age global temperature should reduce enough to trigger 
an IAE. The record of IAEs on Earth is fully consistent 
with the predicted and observed CRF variation — both in 
period and in phase. Moreover, the fit improves when the 
predicted lag in the IAEs after each crossing is included 
and when the actual asymmetric location of the arms is 
considered. Moreover, a random mechanism to generate 
the IAEs is excluded. Nevertheless, one should bear in 
mind that the weakest link still remains the glaciological 
record with its uncertainties. That is, more research on 
the timing and extent of glaciations is required. 

The last agreement is between the Eon time scale star 
formation activity of the Milky Way and presence or com- 
plete absence of IAEs. Here, a more detailed research on 
the SFR activity would be useful to strengthen (or per- 
haps refute) the long term correlation. 

If the apparent correlation between observed CRF vari- 
ations and climate on Earth is not simply a remakable co- 
incidence, an unavoidable question is what is the physical 
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mechanism behind the CRF/temperature relation? Cur- 
rently, there is no single undisputed mechanism through 
which cosmic rays can affect the climate. There are how- 
ever several observational indications that such a relation 
could exist. For example, Forbush events during which 
the CRF suddenly drops on a time scale of days were 
found to correlate with the amount of "storminess" as 
encapsulated by the vorticity area index f30|| , or a concur- 
rent drop in the cloud cover pl[ . There were also claims 
that the galactic CRF, which his modulated by the so- 
lar cycle and slightly lags behind it, correlates with the 
low altitude cloud cover variations [^4], [52). Clearly, an 
in depth study on the possible climatic effects of cosmic- 
rays is imperative. 
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